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Small  Aperture  Analysis  of  the  Dual  TEM  Cell 

and  an  Investigation  of  Test  Object  Scattering 

in  a  Single  TEM  Cell 

P.  F.  Wilson  and  M.  T.  Ma 

Electromagnetic  Fields  Division 

National  Bureau  of  Standards 

Boulder,  Colorado  80303 

Small  aperture  theory  is  used  to  investigate  the  dual  TEM  cell.  Analyzing  coupling 
through  an  empty  versus  a  loaded  aperture  leads  to  a  model  of  dual  TEM  cell  shielding 
effectiveness  measurements.  Small  obstacle  scattering  yields  results  for  both  the  field 
perturbation  and  the  change  in  a  cell's  transmission  line  characteristics  due  to  the 
presence  of  a  test  object  in  a  TEM  cell.  In  each  case,  theoretical  values  are  compared 
to  experimental  data. 

Key   words:     dual   TEM   cell;   perturbation;   shielding   effectiveness;   small 
aperture/obstacle  theory;  test  object  loading. 

1.   Introduction 

The  transverse  electromagnetic  (TEM)  cell  is  well  established  as  a  device  for  creating  a  known, 
broadband,  isolated  test  field  [1].  TEM  cell  applications  include  the  emissions  and  susceptibility 
testing  of  electronic  equipment  for  electromagnetic  compatibility  (EMC)  purposes  [2-8],  the 
calibration  of  antennas  and  probes  [9,10],  as  well  as  the  study  of  biological  effects  of  radio- 
frequency  (rf)  radiation  [11-13].  The  basic  TEM  cell  design,  as  shown  in  figure  1,  consists  of  a 
section  of  rectangular  coaxial  transmission  line  (RCTL)  coupled  via  tapered  sections  to  ordinary 
50  Q  coaxial  cable.  Excluding  the  gap  regions,  the  central  RCTL  test  section  supports  a  linearly 
polarized  TEM  mode  which  is  similar  to  a  free  space  plane  wave.  Both  the  TEM  cell's  waveguide 
properties  [14-20],  especially  the  appearance  of  the  higher-order  modes,  and  the  manner  in  which  TEM 
cell  results  may  be  related  to  the  free  space  environment  [3-5,  21-24]  have  been  carefully  studied. 

As  the  usage  of  TEM  cells  has  become  increasingly  refined,  interest  has  grown  in  applying  its 
well  defined  characteristics  to  additional  metrology  areas.  One  such  example  is  the  dual  TEM  cell 
concept  which  is  being  proposed  as  a  method  to  measure  the  shielding  effectiveness  of  materials 
subjected  to  rf  radiation  [25-29],  As  the  name  suggests,  the  dual  TEM  cell  consists  of  a  pair  of 
cells  coupled  via  an  aperture  in  the  common  wall.  A  comparison  of  the  power  coupled  through  an 
empty  aperture  to  that  transferred  through  the  aperture  when  covered  with  the  test  material  provides 
a  relative  measure  of  the  material's  shielding  effectiveness.  In  general  such  a  waveguide  coupling 
problem  is  quite  difficult  to  analyze.  Fortunately,  the  TEM  cell  is  basically  a  low  frequency 
device.  Therefore  dual  TEM  cell  coupling  is  well  suited  to  an  investigation  based  on  the  small 
aperture  theory  pioneered  by  Bethe  [30].  Although  the  dual  TEM  cell  is  the  primary  motivation  for 
the  present  study,  the  basic  results  may  be  formulated  in  a  fairly  general  fashion  and  thus  are 


applicable  to  other  TEM  cell  associated  problems  such  as  the  loading  effect  on  a  single  TEM  cell  by 
a  small  object. 

Section  2  summarizes  the  basics  of  small  aperture  theory,  as  well  as  the  dual  small  obstacle 
formulation.  If  the  aperture  or  obstacle  is  electrically  small,  then  the  scattering  effect  is 
essentially  equivalent  to  that  produced  by  an  appropriate  set  of  dipole  moments.  These  dipole 
moments  may  be  used  to  predict  the  scattered  fields,  which  in  the  case  of  the  dual  TEM  cell  give  a 
description  of  the  expected  aperture  coupling.  The  dipole  moments  depend  on  both  the  incident 
(exciting)  fields  and  the  shape,  size,  and  orientation  of  the  aperture  or  obstacle.  These  latter 
effects  are  summarized  in  a  quantity  termed  the  aperture  (obstacle)  polari zabi 1 ity .  As  discussed  in 
section  3,  few  polari zabi 1 ities  may  be  derived  directly.  Therefore,  various  approximations  such  as 
numerical  approaches  or  empirical  values  are  inevitable.  Nonetheless,  in  practice  Bethe's  basic 
technique  has  been  shown  to  yield  excellent  results  [31]. 

In  section  4,  small  aperture  theory  is  applied  specifically  to  the  dual  TEM  cell  and  the 
shielding  effectiveness  procedure.  First,  coupling  via  an  empty  aperture  is  considered,  including 
the  case  of  an  appreciably  thick  aperture.  The  agreement  between  theoretical  results  and  measured 
data  is  found  to  be  quite  good  and  serves  to  establish  the  usefulness  of  the  present  approach. 
Coupling  through  a  material -covered  aperture  (loaded)  is  less  straightforward  to  analyze  due  to  the 
more  complicated  geometry  and  the  difficulty  in  specifying  certain  material  parameters.  The 
particular  loaded  aperture  model  chosen  here  does  yield  curves  which  qualitatively  agree  with 
measured  shielding  effectiveness  data.  More  importantly,  our  analysis  provides  insight  into  the 
basic  coupling  mechanism.  In  doing  so,  it  points  out  which  parameters  are  the  most  critical  to  dual 
TEM  cell  measurements  and  thus  which  factors  need  to  be  controlled  most  carefully  in  order  to  obtain 
meaningful  shielding  effectiveness  data. 

Normal  TEM  cell  applications  involve  introducing  a  test  object  into  the  cell.  To  some  extent 
the  test  object  will  perturb  the  TEM  cell  environment  from  its  well  studied  empty  state.  The  TEM 
mode  will  be  scattered  and  higher-order  modes  excited,  although  the  latter  will  usually  be  quickly 
attenuated.  If  the  TEM  cell  is  viewed  as  a  transmission  line  circuit,  then  the  test  object 
represents  a  load  which  will  alter  the  cell's  overall  impedance  characteristics.  Typically,  the 
test  object  size  is  restricted  to  some  fraction  of  the  pertinent  cell  dimensions  and  any  loading 
effects  are  ignored.  Based  on  small  obstacle  theory,  we  are  able  to  derive  an  equivalent  circuit 
representation  for  the  test  object  loading.  The  resulting  loaded  cell  transmission  line  circuit 
yields  the  expected  impedance  change.  In  particular,  impedance  curves  are  generated  for  a  sequence 
of  centrally  located  spherical  conductors  of  increasing  diameter.  Results  suggest  that  the  change 
in  impedance  is  not  likely  to  be  severe  enough  to  cause  problems  in  a  normal  cell  application.  In 
addition  to  the  above  transmission  line  circuit,  section  5  presents  results  for  the  scattering  of 
the  first  four  RCTL  modes.  The  modes  themselves  are  discussed  in  Appendix  A.  It  appears  that  the 
field  variation  from  the  empty  cell  case  can  be  significant,  especially  the  backscattered  TEM 
mode.  Together,  the  loaded  transmission  line  circuit  model  and  the  field  perturbation  results, 
enable  one  to  better  judge  whether  a  particular  measurement  procedure  needs  to  account  for  the 
effects  of  test  object  presence. 


A  short  conclusion,  summarizing  results,  and  suggesting  further  TEM  cell  problems  of  interest, 
is   given   in  the  final    section. 

2.       Small   Aperture   (Obstacle)  Theory 

In  free  space  a  source  distribution  (currents,  charges,  etc.)  may  be  represented  by  an 
equivalent  multipole  expansion  which  generates  the  same  fields  outside  of  some  minimal  volume 
containing  the  sources  [32],  If  the  sources  are  small  compared  to  the  free-space  wavelength,  then 
the  initial  multipole  expansion  terms  are  sufficient  to  predict  the  far  fields.  For  example,  this 
is  the  justification  for  replacing  small  current  loops  with  magnetic  dipoles.  The  same  basic 
argument  holds  in  a  confined  region,  such  as  a  waveguide,  although  with  some  modification  [33],  The 
sources  now  excite  waveguide  modes  or  cavity  resonances.  If  only  the  initial  multipole  terms  are 
retained,  then  the  number  of  modes  whose  excitation  may  be  accurately  described  is  limited. 
Typically  though,  the  higher-order  modes  are  not  of  any  great  interest  because  they  are  rapidly 
evanescent  away  from  the  scattering  region.  Thus,  equivalent  to  the  far  field  restriction  in  the 
free  space  case,  results  should  be  valid  except  in  the  immediate  vicinity  of  the  source  region  where 
the  higher-order  modes  may  still    contribute  substantially. 

Small    aperture   theory,    as   well    as    its    small    obstacle   dual,    is    based   on    generating  the   leading 

terms    in    such    a    multipole    expansion.      Most    details    may    be    found    elsewhere    [31,34,35]    and    need    be 

only    briefly    summarized    here    as    applicable   to   the    TEM    cell.      The    notation    is    essentially    that    of 

Collin       [34],  Given      a      discontinuity      in      a      waveguide      excited      by      some      incident      field 

-       -  -  +     -  + 

distribution  E..   H.  ,  the     scattered     fields  E   ",   H       may     be     expanded     in     terms     of     the     normalized 
11  s  s 

-  +     -  + 
waveguide  modes  E      ,   H       according  to 
n  n 

E   ~   =    V  a     E        ,H       =   I  a     H 
s  Lnns  Lnn 

n  n 
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E  "  =   I  b     E  "         H  "  =   Y  b     H" 
s  nn  n  nn 
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where  +  indicates  propagation  away  from  the  discontinuity  in  either  the  forward  (+)  or  backward  (-) 
direction.  If  we  let  the  z-axis  lie  along  the  direction  of  modal  propagation  and  reference  the  zero 
plane  to  the  discontinuity  center  (arbitrarily  defined),  then  the  normalized  waveguide  modes  may  be 
written  as  follows: 


-  +   ,-     -     ,  +18  z 
E   =  (E   +  a  E   1  e  Pn 
n     nt    z  nz 


H  ±  =    f±  H   +  a  H  )  e+1pnZ 
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where  (3  is  the  propagation  constant  of  the  n-th  mode,  a  time  convention   exp  (ioit)  has  been 

suppressed,  and  the  transverse  field  components   E   and  H   are  related  via  the  admittance 

nt      nt 

dyadic  Y 
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H  .  =  Y   •  E  .  (3) 

nt    n    nt 


Implicit  in  (2)  is  that  the  field  components  inside  the  parentheses,  E  ,  H   etc.,  depend  only  on 
the  transverse  position  (x,  y).  The  orthonormal i zation  condition  takes  the  form 


J  (E   x  H    •  a  ds  =  6  4 

1    y   mt    nt '         z      inn 
CS 

where  CS  is  the  waveguide  cross  section  (assumed  independent  of  z)  and  6   is  the  Kronecker  delta 

mn 

function . 

The  leading  multipole  expansion  terms  in  kQ,  the  free  space  wave  number,  are  the 
electric  (P)  and  magnetic  (M)  dipole  moments.  Given  this  equivalent  source  representation  of  the 
discontinuity,  the  expansion  coefficients  an  and  bp  are  found  to  be  [34] 
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n   2  *■  o  n        n     ; 
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b  =  ■=-"  t\x     H   •  M  -  E   •  P) 
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where   u     is    the    free    space    permeability.       The    dipoles    are    located    at    the    defined    center    of    the 

discontinuity    region.      The    incident    field,    which    excites    the   sources,    is    assumed    to    be   essentially 

uniform    over    the     electrically     small     scattering    discontinuity.        Therefore     it     is     reasonable    to 

decompose    the    dipole    moments    into    the    incident    field    at    the    dipole    location  E    ,   H     times    a    term 

which    characterizes    the    scatterers    geometric   and   material    properties.      These    latter   quantities    are 

the     electric       (a   )   and     magnetic    (a   )  polarizabi 1 ities     and     are     in     general     dyadics.         In     terms 

of   a     and  a     the  expansion  coefficients   become 
em 
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where  the  free  space  permittivity  e  has  been  introduced  for  symmetry. 


This  is  the  basic  equation  describing  waveguide  scattering.   In  order  to  apply  (6)  to  a 
particular  problem,   three  quantities   need  to  be  specified;   (1)  the  normalized  waveguide 

+ 


-  + 
modes  E 
n 


H 


2 


the    polarizabilities     a  and  a  ,  and    (3)    the    incident    field 

e     m 


distribution  E.,  H.  .  The  normalized  waveguide  modes  in  a  TEM  cell  can  be  found  using  either  a 
conformal  mapping  solution  [21]  or  a  singular  integral  equation  approach  [16,23].  Polarizabilities 
constitute  a  problem  of  longstanding  interest  with  solutions  ranging  from  the  relatively  simple 
static  solution  of  Bethe  [30]  to  more  recent  computer  based  numerical  techniques  [36],  Finally,  the 
incident  fields  must  be  detailed.  In  the  case  of  a  dual  TEM  cell  (see  fig.  2),  the  incident  field 
is  the  TEM  mode  in  the  driving  cell  while  the  scattered  fields  may  be  analyzed  in  either  the  sensing 
cell  or  the  driving  cell  itself.  In  discussing  test  object  loading,  the  TEM  mode  will  again  excite 
the  scattering  obstacle  and  the  scattered  fields  are  confined  to  that  cell. 

Before   proceeding   we   may   take   advantage   of   relationships   between   the   modal 

components  E  "  and  H  "  to  expand  (6)  somewhat.  The  admittance  dyadic  Y  is  given  by 
n      n        .  n 
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where  Zn  is  the  wave  impedance  of  the  n-th  mode.   Substituting  (7)  into  (6)  and  assuming  that  the 
polarization  dyadics  are  diagonal  yields 
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This  may  not  appear  to  be  a  useful  expansion  of  (6).  However,  consider  the  following;  if  the 
incident  field  is  a  forward  traveling  TEM  mode  (n=0),  as  will  be  the  case  for  the  two  present 
applications,  then 

H.  -   Y   •  E. 
it    o    it 


H.  =  E.  =  o 

IZ     1Z 


Substituting  these  restrictions  into  (8)  yields 
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where  n  is  the  free  space  wave  impedance.  In  terms  of  the  fields,  we  need  specify  only  in  (10)  the 
electric  field  components  of  the  n-th  mode  and  the  incident  TEM  mode,  at  the  dipole  location.  This 
simpler  form  (10)  will  be  the  basic  equation  needed  for  the  present  investigation.  We  next  consider 
the  polarizabil ities   needed  to  apply  the  above   result. 

3.       Polarizabil ities 

The  electric  and  magnetic  dipole  moments  may  be  defined  in  terms  of  the  currents  and  charges 
induced  on  the  scattering  object  by  the  incident  field  distribution.  If  the  scatterer  is  a  perfect 
conductor  (e  ■>  »)  or  a  perfect  magnetic  conductor  (u.  ■*■<=),  then  P  and  M  may  be  defined  as 
follows   [34] 

(electric   conductor)  (magnetic  conductor) 

~£o 
P  =  /  r   p     ds  or  — ^—     f  r  x  J     ds 

s       e  2      S 


:n: 


M  =         Tr/rxJds  or  —     frpds 

2      }c  e  u.       i        Km 
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where  p^  and  p  are  the  electric  and  magnetic  surface  charges,  J  and  J  are     the  electric  and 
e     pm  3  3     e     m 

magnetic  surface  current  densities,   r  =  xa  +  ya  +  za  is  a  position  vector,  and  S  denotes  the 

x    y     z 

scattering  surface.   For  nonperfect  materials,  that  is   e  and  u  finite,  an  additional  volume 

integral  is  required  [33]. 

Notice  that  for  a  perfectly  conducting  obstacle,  (11)  implies  that  the  electric  dipole  moment 
is  due  to  surface  charge  while  the  magnetic  dipole  moment  is  due  to  circulating  currents.  Thus 
locally  over  a  small  surface  element  ds  the  electric  dipole  moment  is  tangential  to  the  surface 
while  the  magnetic  dipole  is  in  the  direction  of  the  surface  normal.  For  a  perfect  magnetic 
conductor,  duality  suggests  that  these  roles  are    reversed.   In  fact,  if  we  consider  the  fields 

(E  ,  H  )  for    the    electric    source  fp  ,  J  )   case    versus    the   magnetic    source    case 
e   e  v  e   e 

(E   H  ,  P  ,  J  ,)  then    duality    principles    [37]     (E  -  H  ,  H  -  -E  ,  p  -  p  ,  J  -  J  ,  and 
m   m   m   m  J         K  *•  e    m   e    m  ^e   m   e   m 

^o  **  eo^  ^mp^  the  following  relationships 


P  (electric  conductor)  =  e  a  (electric  conductor)  •  E 

o  e  e 

=  J  r   Pe  ds  -  J  r   Pm  ds  (12) 

=   \x     a  (magnetic  conductor)    •  H     =   u.    M   (magnetic  conductor) 
o     m  mo 

and 

M  (electric  conductor)  =  a  (electric  conductor)  •  H 

m  e 

-\    J"rx5eds  .1     /r  x  5mds  (13)  (13) 

=  a  (magnetic  conductor)  •  -E  =  —  P  (magnetic  conductor) 
e  me 

o 


Together  with  the  comments  above,  equations  (12)  and  (13)  yield 


a  (electric  conductor)  -  a  (magnetic  conductor] 
e  m    3 


a  (electric  conductor)  -  a  (magnetic  conductor) 
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In  summary,  the  polarizabil ities  are  not  independent.  Rather,  only  certain  canonical  problems  need 
be  solved. 

It  remains  to  relate  the  above  results  to  a  small  aperture  which  is  neither  a  perfect  electric 
or  magnetic  conductor.  We  find  that  the  polarizabi lity  for  the  above  idealized  aperture  is  related 
to  that  of  a  magnetic  conducting  disk  of  the  same  shape.  In  fact,  the  relationship  is  [34], 


a  (aperture)  =  —  a  (magnetic  disk) 
e  4  e    3 

(15) 

a  (aperture)  =  —  a  (magnetic  disk) 

m  4  m 

A  factor  of  1/2  is  introduced  because  the  disk  has  currents  on  two  sides  while  the  aperture  is  one- 
sided and,  because  the  aperture  dipole  moments  radiate  into  two  half-spaces,  the  effective  dipole 
moments  are  reduced  by  an  additional  1/2  multiplier. 


The  various  integrals  needed  to  evaluate  pol arizabil ities  are  difficult  to  analyze  directly. 
Thus,  a  number  of  approximations  are  employed.  As  mentioned,  the  simplest  approach  is  the  static 
model  detailed  by  Bethe  [30]  which  has  antecedents  in  certain  acoustical  and  electrical  problems 
treated  by  Lord  Rayleigh  [38,39]  and  fluid  dynamics  considerations  given  by  Lamb  [40].  In  practice, 
the  infinite  plane  screen  idealization  requires  that  the  aperture  not  be  located  too  close  to  the 
waveguide  walls.  Fortunately,  if  the  aperture  or  obstacle  is  roughly  a  principal  dimension  away 
from  additional  conducting  walls,  then  the  image  terms  do  not  significantly  affect  the 
polarizabi 1 ity   [41]. 

Conn  [42-43]  has  measured  values  for  a  number  of  interesting  shapes  including  rectangular 
apertures,  rosettes,  and  dumbells,  using  the  electrolytic-tank  method.  More  recently,  the  advent  of 
high-speed  computers  has  led  to  numerical  solutions  based  on  static  assumptions  [44,36],  Rayleigh 
series  (kn)  corrections  to  the  static  solutions  [33,36,45-47],  and  time  domain  considerations 
[48,49],  The  latter  paper  by  Butler  et  al  .  contains  a  very  useful  survey  of  aperture  related 
problems . 

An  alternate  approach  to  computing  the  polarizabil ity  directly  is  to  place  bounds  on  its 
possible  values  [50,51].  Jaggard  and  Papas  [51]  note  that  a  circle  contains  the  maximum  area  for  a 
given  perimeter  length.  Thus,  the  polarizabi 1 ity  of  an  arbitrary  convex  aperture  is  bounded  below 
by  the  value  for  a  circle  with  the  same  area  (minimum  perimeter)  and  above  by  the  value  for  a  circle 
with  the  same  perimeter  length  (maximum  area).  These  bounds  would  be  especially  useful  for  shapes 
of  an  essentially  circular  nature  such  as  an  octagon. 

TEM    cell     applications    will     most    likely    involve    spherical,    or    box-shaped    scattering    objects. 

Apertures   will    also   probably   be   circular,   or   rectangular.     Therefore  only  a   few  of  the  above   results 

need     to     be     explicitly     detailed     here.         The     polarizabi 1 ity     of     an     arbitrary     box     (rectangular 

paral lelopiped)    is    not    readily    available.      However,    the    above    logic    suggests    that    bounds    based   on 

spheres    of    equal     surface    area    (above)    and    volume    (below)    should    give    an    indication    of    the    proper 

result.         For     the     types     of     scatterers      considered     here     only     the     diagonal      elements     of     the 

polarizabi 1 ity    dyadics    are    nonzero.       Letting    the    three    orthogonal    axes    have    unit    vectors      x.,  we 

need  only  specify     a    .   =   x.    •   a      •  x.   and      a    .   =   x.    •   a      •  x.. 
ejjej  mjjmj 

For  ellipsoidal    conducting  obstacles  we   find  that  the  static  approximation   gives    [34] 

a    .   =   V/L. 
ej  J 

(ellipsoid)  (16) 

a    .    =   V/(L.    -    1) 
mj  J 

4 
where  V   =  T7triroro      is    the    volume,    r-,  ,    r0,    and    r-j    are    the    semi-axes,    a   .   and      a   .   are    associated 
3123  '      1 '      c '  3  ej  mj 

with  the   respective  semi-axes,   and  the   integral    L-   (j   =   1,   2,    3)    is  defined   by 


J  2  j  „  0  0  0        1/2  {U) 


0     (s  +  r/)[(s  +  r^)(s   +  0(s  +  r/)    ] 


The  L-  integrals  may  be  evaluated  in  terms  of  elliptic  integrals.   In  the  case  of  a  sphere  (r-,  =  r-o 
=  r^  =   r)    L-  =  1/3,  and  the  pol ari zabi 1 i ties  reduce  to 

A       3 
a  .  =  4-nr 

ej  (sphere)   .  (18) 

a  .    o  3 
mj  =  -2itr 

If   we    let    rj    =    r^   and    r^   =    0,    then    the    result    is    a    circular  disk   of    radius    r.      Equations    (16)    and 

(17)   then  yield     a  ,    =  ao  =  ao  =  0,    and 

J  ml  m2         e3 

16     3 
ael=aeZ=Tr 

_o     3  (circular  disk)      .  v      ' 

a  -.  =  — ~  r 
m3         3 

In  the  case  of  a  square  conducting  disk  of  side  length  d,  the  polarizabil  ities  are  approximately 
given  by  [52] 

3 
a  , =  a  0  =  1 .032  d 
el   e2 


a  -  =  -0.455  d 
m3 


i  ,-n  (20) 

o        (square  disk) 


The  polarizabil  ities   a  , ,  and   a  ?  are  associated  with  axis  aligned  with  the  sides  of  the  square 

while   a  -  is  associated  with  the  surface  normal  direction.    As  with  the  circular  disk, 
eo 

ml    m2    e3 

Using  duality  (15),  equations  (19)  and  (20)  immediately  give  the  pol arizabi 1 ity  for  a  circular 

and  square  aperture.  Specifically  we  find  that  a  ,  =  a  „  =  a  ~  =  0, 


2     3 
a  0  =   — o  r 
e3         3 


a  ,  =   a  o   =  sr  r 

ml       m2       3 

and 

a  -   =   -0.114d3 
e3 


.     _  (circular  aperture) 


^  (square  aperture)      .  (22) 


a   ,  =  a  0  =  0.258d 
ml    mi 

In  the  future,  we  will  simply  speak  of  the  aperture  polarizabilities   a  and   a  .    Their 
orientations  (normal  and  tangential)  will  be  understood. 

The  validity  of  (21)  and  (22)  is  based  on  the  assumption  of  an  infinitely  thin  conducting 
plane.  In  fact,  the  aperture  could  have  an  appreciable  thickness.  Cohn  [53]  and  others  have 
suggested  that  a  thick  aperture  be  modeled  as  a  section  of  waveguide  operating  below  cutoff.  The 
equivalent  dipole  moments  effectively  excite  the  evanescent  modes  in  the  aperture  waveguide.  These 
modes  then  generate  a  new  set  of  dipole  moments  at  the  aperture's  other  end.  A  very  simple 
thickness  correction  follows  from  only  considering  the  lowest-order  mode  excited  individually 
by  P  and  M  .  For  example,  given  a  square  aperture,  the  electric  dipole  moment  P  is  normal  to  the 
aperture  and  thus  longitudinal  to  the  aperture  waveguide.   Therefore,  P  will  excite  TM  modes  with 
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the  TM-jj   being  the   first  mode  for  a   square  cross   section.     Thus    a     should   be   replaced   by 


TM 

a    e"PllT  (23) 

e 


TM 
where  3  '  is  the  propagation  constant  of  the  TM-q  mode  and  T  is  the  aperture  thickness, 

Likewise,  M  will  excite  both  TE  and  TM  modes,  thus  the  TE1(j  mode  effect  will  be  included, 

Therefore  a     is  replaced  by 
m     K      J 

-pJot 

(24) 


where  pin  is  the  TEiq  propagation  constant.  Alternate  thickness  corrections  are  provided  by 
McDonald  [54]  who  uses  the  reaction  concept  to  generate  an  equivalent  circuit  for  the  thick  aperture 
and  by  Akhiezer  [55]  who  generalizes  Bethe's  approach. 

Shielding  effectiveness  measurements  involve  loading  the  aperture  with  the  material  under 
test.  If  the  aperture  is  actually  fitted  with  the  test  material,  then  the  above  waveguide  model  can 
be  applied.  One  simply  modifies  the  propagation  constants  to  reflect  the  new  waveguide  medium. 
Realistically,  it  will  be  more  convenient  to  cover,  rather  than  plug  the  aperture  [25-26].  Thus,  an 
alternate  approach  to  determine  the  loaded  polarizabi 1 ities ,  designated  a  and  a  ,  is  needed. 

Some  authors  have  generalized  the  basic  aperture  problem  to  allow  for  different  media  on  each 
side  of  the  conducting  plane  [56-58],  The  static  approximation  yields  [56]. 

2e, 


a  /a 


e      e        e,    +    e9 


a  /a     = 


(different  media)  (25) 


mm       u,    +   \u 


where  j  =  1,  2  determines  from  which  side  the  incident  field  originated  and  the  direction  of  the 
electric  dipole  moment  is  into  the  material  on  each  side  of  the  aperture.  These  values  might  be 
useful  if  the  sample  covering  was  very  thick.  However,  in  practice,  thin  samples  are  more  likely. 

Therefore  the  problem  of  a  thin  sheet,  possibly  lossy,  is  of  interest.  Marin  [59]  gives  some 
static  approximations  for  the  polarizabi 1 ity  of  a  circular  aperture  when  a  dielectric 
sheet  (a  =  0)  of  thickness  h  is  present.  He  derives  closed-form  expressions  for  very  thick  and  very 
thin  sheets.  In  the  latter  case  for  example,  he  finds  that 

~  /     i  ,  3^Er2  "  X)  r2hx  „  r2h>  .   r2h »  ,9,* 

ae/ae  =  1  +       2TO9         (— )   *"(— )  +  °(— )  (26) 

r2 

(dielectric  sheet-circular  aperture   radius   r) 
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where  e  „  is  the  dielectric  constant  of  the  material  and  2h/r  is  the  quantity  assumed  small  in  the 
thin  sheet  case.  An  air  backing  is  assumed  on  each  side  of  the  sheet.  As  may  be  seen,  the 
normalized  pol arizabil ity  (26)  goes  to  unity  as  the  thickness  h  vanishes.  The  magnetic 
pol arizabi 1 ity  will  be  unaffected  by  the  dielectric  sheet  and  therefore  is  still  given  by  the  second 
expression  in  (21) . 

Actual  shielding  materials  will  be  conductive  (a  »  1)  since  it  is  desirable  that  a  shield 
material  both  reflect  and  absorb  fields  as  well  as  provide  a  continuous  path  for  induced  surface 
currents  [60].  A  lossy  sheet  presents  a  further  complication  in  determining  the  pol  arizabi 1 ities . 
In  the  static  limit,  the  sheet  will  appear  as  a  perfect  conductor  to  the  electric  fields  and 
transparent  to  the  magnetic  fields  [61],  that  is  a  ■*  o  while  a  remains  unchanged.  These  zero-th 
order  terms  provide  some  insight  into  the  expected  behavior  of  a  loaded  aperture;  however,  the 
dependence  on  the  material  properties  is  still  hidden. 

Casey  [60]  has  developed  an  integral  equation  solution  for  the  case  of  wire  mesh  loading.  By 
modeling  the  mesh  in  terms  of  an  equivalent  resistive  sheet,  he  is  able  to  generate  a  first  order 
correction  to  the  static  limit  values.  As  with  Marin's  effort  [59],  Casey's  results  are  for  a 
circular  aperture.  He  finds  the  following 

- 

a  /a   =  [1  +-— ] 
e  e       3u 


r-,      In 

a  /a  =  [1  +  ,  ..    ~ 

m  m       3it(l  +  z 
v     c 


'wire  mesh  -  circular  aperture  radius  r)  [21] 

,-1 


where 


= 

1  um   r  • 
^0 

m 

21 

s 

C     r 
s 

e 

2e     a 
o     s 

zc  =   2*yZs  (28) 

Zs   =   Zw  as  +   iwLs 


-2ur  /a     -1 
in   (1  -  e         u    s) 


Cs        2«e0(erl   +   er2)    a$ 


n  a  -2-nr  /a     -1 


Ls'-^O-e       "   S) 
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ac  denotes  the  mesh  size,   r   is  the  mesh  wire  radius,  Rr      is  the  contact  resistance, 

S  CO  c  * 

and  e  ,  and  e  0  are  the  dielectric  constants  of  the  material  (if  any)  on  each  side  of  the  mesh.  The 
internal  impedance  per  unit  length  Z'  is  given  by 

A  I  (A) 
co     oo  2  I,  (A) 

where  R1    is   the  dc    resistance   per   unit    length   of  the  mesh  wire,      A  =    Moo  x    J        ,    x     =   u.  o  r       is  the 
to  ^uroocoojoj 

diffusion    time    constant    of   the   wire,      \x     and      a    are   the   wire's    permeability    and    conductivity,    and 

co        to  r  J 

I  (y  =  0,1)  denotes  a  modified  Bessel  function  of  the  first  kind.  These  expressions  are  somewhat 
formidable  but  are  included  because  many  shield  materials  are  indeed  such  wire  mesh  loaded  sheets. 

The  polarizabil ities  in  (27)  also  apply  to  a  simple  lossy  sheet  if  we  properly  modify  the 
equivalent   sheet   impedance   Z,.   in   (28).     If   we   let    r  /a.  +0  (zero   wire   thickness) 

^  CO   S 

then  C<-  -*  °°  and  L<-  -*  0  implying  that  Zs  has  only  a  resistive  term  and  that  p  ->  =°,  which  means  the 
normalized  electric  polarizabi 1 ity  in  (27),  is  again  zero  as  in  the  static  case.  In  addition,  one 
may  show  that  for  an  electrically  thin  resistive  sheet  [34] 


z   =4t  (3o; 

s   oh  v 


where  a  is  the  conductivity  and  h  the  sheet's  thickness.     Thus  we  have 

i    4u)  li  r   a  h 

am/am    "   t1   +   3n  (1   +   2u  R      ah)^ 

c 

(lossy  sheet  -  circular  aperture  radius  r)  (31' 

a  /a   =0   . 
e  e 


If  the  contact  resistance  R  ,  and  the  material  properties  a   and  h  are   known,  then  these  expressions 
may  be  readily  applied,  as  will  be  seen  in  the  case  of  the  dual  TEM  cell. 

4.   The  Dual  TEM  Cell 

The  basic  goal  of  EMC/EMI  engineering  is  to  insure  that  electronic  equipment  will  function 
properly  in  a  specified  electromagnetic  environment.  This  requires  both  that  a  system  not  be  overly 
susceptible  to  the  expected  ambient  fields  and  that  the  system  itself  does  not  contribute  undue 
amounts  of  spurious  radiation.  Satisfying  these  requirements  has  led  to  the  general  concept  of 
shielding,  as  well  as  numerous  interference  control  guidelines,  test  facilities,  and  test  procedures 
[62-65]. 
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Roughly  speaking,  EMC  test  objects  may  be  separated  into  three  levels;  (1)  materials,  (2) 
components,  devices,  etc.,  and  (3)  systems.  The  desire  at  each  level  may  be  to  study  shielding 
effectiveness  (SE),  however,  each  grouping  presents  a  quite  different  test  problem.  For  example, 
the  SE  of  a  plain  sheet  of  material  may  be  quite  different  from  the  SE  of  a  container,  or  box  formed 
from  the  same  material  [66,67].  Similarly,  a  complex  system  may  fail  to  satisfy  a  set  of  EMC 
requirements,  even  though  each  component  individually  meets  those  same  requirements  [68].  Thus 
although  a  knowledge  of  a  material's  SE  is  of  great  value  in  predicting  its  behavior  in  a  larger 
system,  SE  figures  should  not  be  used  in  an  absolute  sense.  In  addition,  there  is  great  variance 
among  test  results  just  at  the  materials  level  [29]. 

For  materials,  a  SE  test  procedure  based  on  the  dual  TEM  cell  has  certain  advantages  over  the 
widely  used  shielded  room  technique  based  on  MIL  STD-285  [25,26].  Dual  TEM  cell  measurements 
require  less  time,  offer  excellent  repeatability,  are  readily  standardized,  and  require  only 
moderate  sample  sizes.  Perhaps  more  importantly,  given  the  wide  variation  in  test  results  mentioned 
above,  the  dual  TEM  cell  is  more  readily  analyzed  theoretically  than  most  techniques.  This  greatly 
enhances  the  interpretation  of  SE  test  results  since  it  is  clearer  what  is  being  measured  and  what 
values  should  be  expected  for  a  given  material  sample.  Dual  TEM  cell  disadvantages  are  primarily 
the  limited  frequency  range  (low)  inherent  to  TEM  cells  in  general  and  the  difficulty  in  controlling 
the  aperture  environment  and  the  mounting  of  the  test  material.  This  latter  difficulty  is  shared  by 
most  methods  using  limited  sample  sizes  [29],  including  the  new  ASTM  proposed  transmission  line  and 
dual-chamber  standards  [69]. 

As  shown  in  figure  2,  a  dual  TEM  cell  consists  of  a  pair  of  TEM  cells  coupled  via  an 
aperture.  The  aperture  serves  to  transfer  power  from  the  driving  cell,  shown  here  as  the  upper  cell 
with  the  signal  source  located  at  port  1,  to  the  receiving  or  sensing  cell,  here  the  lower  cell  with 
the  output  measured  at  either  port  2  or  port  3.  The  insertion  loss  (IL)  provided  by  the  shielding 
material  is  defined  as  the  ratio  of  the  received  power  (P)  for  an  empty  aperture  to  that  (P),  for 
the  loaded  aperture  case,  or 

IL  =  10  log  |4|  (dB)  .  (32) 

P 

In  each  case  the  input  power  should  be  held  constant.  Received  power  may  be  measured  either  at  port 
2  (backward  port)  or  port  3  (forward  port)  with  the  other  terminated  in  a  50  Q  load,  as  shown  in 
figure  2  for  a  forward  port  measurement. 

Dual  TEM  cell  coupling  may  be  described  via  the  small  aperture  theory  developed  in  the  previous 
two  sections.  Because  a  TEM  mode  is  used  as  the  incident  field,  equation  (10)  applies,  with  (EQX, 
E0  )  being  the  TEM  mode  electric  field  components  in  the  driving  cell  and  (Enx,  E  )  denoting  the 
n-th  mode  electric  field  components  in  either  the  driving  or  receiving  cell.  Although  the  cells  may 
be  kept  general,  we  will  consider  the  simpler  case  of  identical  cells.  In  addition,  if  the  aperture 
is  longitudinally  centered,  then  the  evanescent  modes  should  not  couple  appreciable  power  to  the 
output  ports.  Thus  only  the  TEM  mode  expansion  coefficients  aQ  and  bQ  are  of  interest.  For 
identical  cells,  equation  (10)  yields 
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'33) 


-i    k 

a      =      -      °    {(a        +   a      )    E        E.       +    (a        +   a      )    E        E.     } 
o         2rj  ex         my       ox     ix  ey         mx       oy     iy 


-i  k 

b   =  -=— -  {(a    -  a   )  E    E.   +  (a   -  a   )  E    E.  } 
o    2ti     ex    my'  ox  ly     ey    mx   oy  ly 


The  normalized  field  components  in  (4)  for  the  TEM  mode  are  found  in  Appendix  A  (A-10). 
Referring  to  figure  3,  which  defines  the  relevant  RCTL  dimensions  (a,  b-,  g)  as  well  as  the  x-y  axes 
orientation,  we  find  that  for  symmetrical  cells  (b^  =  b£  =  b) 


r         2   ,1/2     r       smh  M   (b   -  py)      .      M  .      .,      ,    ,M  , 

E     =  —  Z           )       — u— ~-r — iLLJ-  sin  Mx  sin  Ma  J    ( Mg , 

xac          L,        smhMb  ova' 
mo=l 


!34) 


P  2   ,1/2     r       cosh   M   (b  -  py)  .,  M    . ,    /,.  » 

E     =  p  —  Z  )       •   u  m  k"^0-  cos  Mx  sin  Ma  J    (Mg) 

yKac  L  .  sinn  M  b  o     3 

J  mo=l 


where  the  characteristic   impedance   1Q  is   given  by 

-1 

Z     =   t,    |     [in    (--)    -       X      -  (1   "  "th  Mb)]  .  (35) 

c         o   8  iig  ,   m 

mo=l 

mo  refers  to  the  summation  over  odd  m  only  (m=l,3,5  .  .  .),  M  =  mu/2a,  and  p  is  a  sign  function, 
with  p=l  above  the  septum  and  p=-l  below  the  septum.   In  the  driving  cell  E^x  =  Ex  and  E^   =  E 
while  in  the  receiving  cell  EQX  =  Ex  and  E   =  -E  .   Notice  also  that  (x,  y)  need  to  be  properly 
defined  in  each  cell.  For  example,  the  aperture  is  located  at  y=-b  in  the  upper  cell  but  at  y=b  in 
the  lower  eel  1  . 

Insertion  loss  measurements  are  based  on  the  power  coupled  to  the  receiving  cell.  Thus,  we 
will  concentrate  on  the  expansion  coefficients  a  and  b  in  the  lower  cell.  If  the  aperture  is 
centered  transversely,  (x=0)  .  then  E   =  E.  =  0  and  equations  (33)  and  (34)  combine  to  give 

o  /         \  r  2 
a   =  -s —  (a   +  a  )  E 
o    2ti    ey    mx'  oy 

i  k  „ 

b  =  y—2.  (a   -  o  ■)•  E  •  (36) 

o   2nQ    ey    mx'   oy 

r  2  ,1/2   v   sin  Ma   ,  ,.,  .. 

E  =  —  Z  '    )    .  ,   '■  J  (Mg) 

oy  a  c    L   ,  sinh  Mb  o  3 

J  mo=l 
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This  is  the  basic  equation  describing  dual  TEM  cell  coupling.   As  may  be  seen,  the  result  is  quite 

simple.   Notice  that  the  normalization  relationship  in  (4)  implies  that  the  power  carried  by  the 

,2  ,? 

scattered  TEM  mode  is  given  by   a  \'/n   and    b   /0  in  the  forward  and  backward  directions, 

■  o '   2      '  o '   2 

respectively.  Finally,  the  differing  signs  in  the  equations  for  a  and  b  imply  that  power  is 
coupled  asymmetrically  to  the  sensing  cell.  Because  a  and  a  are  of  opposite  sign  for  an  aperture, 
bQ  will  be  greater  in  magnitude  than  aQ  and  thus  the  dual  TEM  cell  acts  more  as  a  backward  than 
forward  coupler. 

The  dual  TEM  cell  used  for  measurements  at  NBS  featured  a  square  aperture,  thus  the 
polari zabi 1 i ties  for  the  unloaded  case  are  given  by  (22).  Substituting  these  into  (36)  yields 

i  k 

a„  =  -5-2-  (0.144)  d3  Ert2 
0    2r|  oy 


'37' 


b  =  -T—2-   (0.372)  d3  E  2 
o    2-n  ■  v     '  oy 


The  relevant  cell  dimensions  are:  a  =  9  cm,  b/a  =  2/3,  g/a  =  0.24,  and  d/a  =  0.56  which  gives  each 
cell  a  characteristic  impedance  Zc  of  approximately  50.6  Q.  Equation  (37)  implies  that  the 
asymmetry  in  coupling  should  be  about  8.3  dB .  We  will  first  employ  (37)  to  predict  dual  TEM  cell 
coupling  for  the  unloaded  case  in  order  to  establish  the  usefulness  of  the  small  aperture 
equations.  Next,  the  loaded  polari zabi 1 ities  developed  in  section  3  will  be  used  to  analyze  actual 
IL  measurements  made  on  a  particular  sample  material.  Finally,  dual  TEM  cell  results  will  be 
related  to  those  of  a  shielded  room  when  coplanar  loop  antennas  are  used,  as  was  done  by  the  Tecknit 
group  [25,26]. 

For  an  unloaded  aperture,  the  power  coupled  to  the  receiving  cell  normalized  by  the  input  power 
is  given  by 


P  =  20  log  -77   la  I     (forward  -  port  3) 

3  2  '  o ' 

or  (381 

P  =  20  log  4  lb  |     (backward  -  port  2) 
a  2  '  o '  K 


Figure  4  shows  the  theoretical  values  for  both  ports  versus  the  measured  values  for  port  3.  Data 
for  port  2  were  not  available  for  this  particular  measurement  run.  The  cells  were  screwed  tightly 
together  assuring  good  contact  and  a  thin  aperture.  As  may  be  seen,  theoretical  and  measured  data 
agree  to  within  a  decibel  over  the  usable  frequency  range  of  these  cells,  which  is  up  to  about  833 
MHz,  where  the  first  significant  higher-order  mode  (TE-^g)  appears.  At  800  MHz,  kQd  =  .84,  therefore 
the  aperture  is  no  longer  electrically  very  small.  However,  results  remain  surprisingly  good 
largely  because  the  TEM  mode  remains  uniform  over  the  aperture. 
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A  second  set  of  measurements  is  shown  in  figure  5  against  the  same  theoretical  curves.  In  this 
case  power  was  measured  at  both  sensing  cell  ports,  again  for  an  unloaded  aperture,  but  the  cells 
were  no  longer  screwed  together.  Instead,  conductive  gasketing  was  placed  around  the  aperture  and  a 
50-pound  weight  was  used  to  create  electrical  contact.  This  configuration  has  the  advantage  of 
quick  access  to  the  aperture  if  materials  are  to  be  tested,  but  does  not  provide  quite  as  clean  an 
aperture  environment  as  when  the  cells  are  screwed  together.  Figure  5  indicates  that  agreement  is 
reasonable,  to  within  a  few  decibels,  but  not  as  good  as  with  the  tightly  joined  cells.  Because  the 
gasketing  significantly  thickened  the  aperture  to  about  6mm,  (23)  and  (24)  may  be  used  to  modify  the 
polarizabi 1 i ties .  The  results  are  shown  in  figure  6.  Agreement  has  been  improved,  especially  for 
the  port  3  measurements.  Taken  together,  these  curves  indicate  the  importance  of  carefully  joining 
the  cells  and  point  toward  the  aperture  environment  and  sample  mounting  as  potential  problems  that 
must  be  carefully  addressed. 

Insertion  loss  measurements  involve  a  loaded  aperture.  The  appropriate  polarizabi 1 ities  for  a 
resistive  sheet  of  conductivity  a,  thickness  h,  and  contact  resistance  Rc  are  given  by  (31)  for  a 
circular  aperture.  Unfortunately,  the  equivalent  result  for  a  square  aperture  is  not  available. 
However,  we  may  generate  an  approximate  loaded  square  aperture  pol  ari zabi 1 ity  by  defining  a 
radius  r  which  yields  the  same  magnetic  polari zabi 1 i ty  a  for  the  unloaded  case,  that  is 


3   4  ~3 
a  =  0.258  d  =  ■=■  r   (equivalent  circular  aperture] 
m  3      M  r 


:391 


Solving,  we  find  that  r  ~  .579  d.  Because  this  equivalent  radius  reproduces  the  unloaded  square 
aperture  polarizabi 1 ity ,  it  should  yield  reasonable  results  for  the  loaded  polarizabil i ty . 
Recalling  the  definition  of  insertion  loss  (32),  we  find  the  following 


IL  (P3)  =  20  log 


IL  (P2)  =  20  log  - 


:40! 


where,  as  in  section  3,  a  tilde  is  used  to  designate  loaded  quantities.  If  we  insert  the  unloaded 
aperture  polari zabi 1 ities  (22),  and  the  loaded  values  based  on  (31)  with  the  equivalent  radius 
defined  above  (39),  then  the  IL  equations  become 


IL  (P3)  -  -5.1  +  20  log 


1  + 


i  4io  li  r  oh 
_q 

3  n  (T  +  2  n  R  IF 
v        c 


(cm; 


IL  (P2)  -  3.2  +  20  log 


1  + 


i  4u>  li  r  oh 

M3 

3n  (T+  2it  R   oh) 
c 


(41) 


:<jb) 
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Note  that  (41)  implies  that  as  the  frequency  decreases  (to  ■>  o)  the  insertion  loss  seen  by  the 

forward  port  (P3)  may  be  negative,  or  more  power  is  coupled  into  the  sensing  cell  when  the  material 

is  present  than  for  the  unloaded  aperture.   The  cause  may  be  seen  in  the  basic  coupling  equation 

(36).    In  the  forward  direction  the  electric  and  magnetic  fields  interfere  destructively 

(a   +  a  ) .  If  a  lossy  test  material   is  introduced,  the  electric  field  is  screened  out, 
„  ey    mx  „ 

a  ■*  o.  Thus  although  a  is  also  reduced,  the  lack  of  interference  may  more  than  compensate  for 
reduced  magnetic  coupling,  and  thus  the  apparent  gain  due  to  the  material  loading.  This  same  effect 
has  been  observed  in  shielded  room  IL  measurements  as  well  [70]. 

Tests  were  performed  on  a  piece  of  sample  material  which  consisted  of  three  layers;  polyester, 
aluminum,  and  vinyl.  The  layers  were  of  thickness  0.01,  0.01,  and  300  micrometers  respectively. 
The  gasketed  aperture  of  figures  4  and  5  was  used  to  seat  the  material.   In  order  to  apply  (41),  we 

need  to  specify  o,  h,  and  R  .  The  dielectric  layers,  polyester  and  vinyl,  will  not  provide 

-8  7 

significant   shielding,   thus  only  the  aluminum  layer   (h  *  10"  m,  o  «  3.7  x  10  S/m)  will   be 

considered.  The  contact  resistance  Rc  for  our  test  runs  is  unknown.  Thus  curves  for  representative 

values   (R  =  0,  1,  2  Q)  will  be  presented.  Figures  7  and  8  show  the  measured  IL  at  the  forward  and 

backward  ports  respectively.   At  the  forward  port,  the  R  =  1Q  curve  well  fits  the  data  while  at 

the  backward  port,  the  measured  data  falls  between  the  R  =  1Q  and  R  =  2Q  curves.   The  important 

point  is  not  to  try  and  match  one  of  these  curves,  since  R  is  unknown  and  the  contact  impedance  may 

in  fact  be  capacitive  rather  than  resistive;  rather,  that  the  curves  qualitatively  show  the  proper 

behavior.   If  indeed  the  contact  resistance  is  known  (measured),  then  (41)  should  yield  accurate  IL 

predictions.   Also  note  in  figure  7  that  as  the  frequency  approaches  a  few  megahertz,  the  forward 

port  IL  does  indeed  go  negative  as  predicted.   In  fact,  in  both  cases,  the  measured  data  tend  toward 

the  proper  static  limit  (41),  that  is  IL(P3)  +-5.1  dB  and   IL(P2)  ->3.2  dB .  Because  a  thick 

aperture  was  used  we  may  again  introduce  the  waveguide  corrections  (23)  and  (24).   The  thickness 

corrected  results  are  shown  in  figures  9  and  10.   The  correction  to  the  previous  figures  is  only 

slight,  since  the  IL  ratios  (41)  tend  to  divide  out  the  thickness  modifications.   Thus  IL 

measurements  should  be  fairly  insensitive  to  aperture  thickness  (a  definite  advantage),  and  this 

effect  may  be  ignored. 

Finally,  we  may  make  some  comments  about  the  comparison  between  SE  measurements  made  in  a 
shielded  room  versus  a  dual  TEM  cell.  If  a  shielded  room  test  procedure  uses  coplanar  loops  to 
transmit  and  receive  via  the  coupling  aperture  [25,26],  then  we  expect  a  low-impedance  field  at  the 
aperture  for  low  frequencies,  that  is  a  largely  magnetic  field.  The  loops  act  essentially  as 
magnetic  dipoles.  Thus,  coupling  is  essentially  through  the  magnetic  pol arizabi 1 ity  alone.  As  in 
(40),  the  geometry  and  antenna  considerations  basically  divide  out  when  determining  IL,  therefore  we 
should  approximately  have 


IL  (shielded  room)  ~   20  log 


(42: 


for  the  case  of  coplanar  loops.  For  other  types  of  antennas,  the  fields  exciting  the  aperture  will 
be  different  and  consequently  so  will  the  polarizabil ity  combination.  If  the  shielded  room  IL  (42) 
is  compared  to  the  forward  port  dual  TEM  cell  IL  (40),  we  find  that 
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IL    (shielded   room)   -    IL   (P3)    -  5.1   dB 


'43: 


The  Tecknit  group  performed  precisely  this  comparison  [25,26],  and,  in  general,  the  shielded  room 
gives  a  greater  IL  as  predicted  by  (43). 

Ondrejka  and  Adams  [29]  present  IL  measurements  on  the  above  described  sample  for  five  types  of 
SE  test  procedures;  (1)  time  domain  above  a  ground  screen,  (2)  time  domain  through  a  loaded  aperture 
in  a  screen  room,  (3)  cw  through  a  loaded  aperture  in  a  screen  room,  (4)  dual  TEM  cell,  and  (5) 
coaxial  transmission  line  with  flange  holder.  The  results  are  shown  in  figure  11,  with  the  curves 
labeled  according  to  the  above  numerical  listing  and  IL  is  defined  as  in  the  dual  TEM  cell  case 
(32).  Clearly,  there  is  significant  variation  among  test  results  and  a  number  of  contributing 
factors  are  discussed  in  their  paper.  However,  the  dual  TEM  cell  curve  (#4)  gives  significantly 
less  IL  than  the  other  procedures  which  are  roughly  grouped  into  a  15  dB  band.  Curve  #4  was 
generated  with  the  thick  aperture  configuration,  consisting  of  the  same  data  as  presented  in  figure 
10  (backward  port).  What  figure  11  demonstrates  is  the  importance  of  controlling  the  contact 
resistance.  If  we  place  the  theoretical  zero  contact  resistance  curve  from  figure  10  onto  figure 
11,  as  shown  by  curve  #6,  then  the  predicted  dual  TEM  cell  IL  is  more  in  agreement  with  the  results 
from  other  techniques,  especially  at  the  higher  frequencies.  The  contact  resistance  is  very 
important  in  dual  TEM  cell  measurements  since  there  are  direct  current  paths  between  the  input 
(transmitter)  and  output  (receiver)  ports.  In  the  other  techniques,  except  for  the  coax  line,  the 
receiver  and  transmitter  are  electrically  isolated,  therefore  the  material  mounting  may  not  be  as 
critical.  For  the  coax  line,  the  large  flanges  which  hold  the  material  create  capacitive  coupling 
which  may  be  sufficient  to  essentially  short  across  the  material  thereby  effectively  giving  a  zero 
contact  resistance.  Clearly,  the  aperture  environment  is  a  critical  factor  which  deserves  continued 
attention . 

The  small  aperture  analysis  of  dual  TEM  cell  coupling  presented  in  this  section  appears  to  be  a 
useful  technique  for  analyzing  SE  and  IL  measurements.  As  the  dual  TEM  cell  concept  is  refined,  the 
above  analysis  may  well  allow  one  to  accurately  predict  the  expected  IL  given  the  material 
parameters  a,  h,  and  R  .  This  may  well  aid  in  theoretically  designing  shield  materials,  wire  meshes 
etc.,  for  specific  shielding  problems.  Alternately,  if  a  given  material's  electrical  properties  are 
unknown,  then  dual  TEM  cell  measurements  could  be  used  to  generate  an  equivalent  sheet  impedance 
which  would  help  in  predicting  the  material's  expected  behavior  in  other  EMI  environments. 

One  final  comment:   insertion  loss  measurements  made  at  both  ports  may  be  used  in  (40)  to  solve 

for  a   and  a  .  This  allows  us  to  consider  the  "magnetic  field"  insertion  loss  (ILm) 
mx      ey  3  m 


IL  =  20  log 
m       3 


mx 


mx 


and  the   "electric  field"   insertion   loss    (ILe) 
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IL  =  20  log 
e       3 


ey 


ey 


separately.   This  should  aid  in  relating  the  dual  TEM  cell  SE  measurements  to  results  obtained  via 

other  methods  [71].    We  have  already  mentioned  that  small  loops  (magnetic  dipoles)  should 

essentially  give  Il_m  (42)  in  the  near  field.   For  conductive  shield  materials,  we  should  have 

IL  <  IL  ,  that  is,  less  shielding  of  the  magnetic  fields.   The  various  SE  test  procedures  should 
me  K 

give  insertion  loss  values  somewhere  between  these  two  limits.  Thus  dual  TEM  cell  measurements  can 
be  used  to  determine  a  shielding  range  expected  for  a  material. 

5.   Small  Object  Scattering  in  a  TEM  Cell 

As  indicated  in  the  introduction,  most  TEM  cell  measurements  involve  placing  some  test  object 
inside  the  cell  and  then  exciting  the  TEM  mode.  If  the  test  object  is  kept  small  compared  to  the 
test  chamber  height  and  its  dimensions  are  much  less  than  the  operating  wavelength,  then  it  is 
reasonable  to  assume  that  the  cell  environment  is  not  significantly  perturbed  from  its  empty 
state.  The  small  aperture  theory  developed  in  section  2  may  be  used  to  investigate  the  effects  of 
test  objects  from  two  viewpoints.  First,  when  only  the  scattered  TEM  mode  is  considered,  the  TEM 
cell  may  be  viewed  as  a  transmission  line  circuit  with  the  test  object  introducing  some  type  of 
load.  Analyzing  the  reflected  and  transmitted  TEM  mode  enables  us  to  represent  the  test  object 
loading  in  terms  of  an  equivalent  circuit.  Thus  the  change  in  the  overall  TEM  cell  transmission 
characteristics  due  to  the  test  object  presence  can  be  studied.  Second,  the  excitation  of  the 
initial  higher-order  modes  may  be  modeled.  This  enables  us  to  theoretically  assess  the  expected 
field  perturbation  due  to  test  object  scattering.  In  each  case,  the  indication  is  that,  indeed,  the 
TEM  cell  environment  is  not  greatly  disturbed  for  the  simple  types  of  scattering  objects  considered 
here.  If  the  sensitivity  of  the  desired  measurement  is  such  that  even  small  deviations  in  cell 
characteristics  are  important,  then  the  present  theory  should  be  useful  for  predicting  what  limits 
should  be  placed  on  the  test  object  size  and  the  measurement  frequencies. 


As  in  the  previous  section,  we  will  again  assume  that  the  scattering  obstacle  is  excited  by  an 
incident  TEM  mode.   Therefore,  the  basic  equation  is  (10).   The  scattered  TEM  mode  (n=0)  will  be 


considered  first. 
(34),  we  find  that 


Noting  that 


E   =  E.  =  E  and  E 
ox    lx    x     oy 


E. 

iy 


E  ,  where  E   and  E   are  given  by 
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The  expansion  coefficient  b  represents  the  reflected  TEM  mode  while  aQ  gives  the  scattered  TEM  mode 
in  the  forward  direction.  Thus,  the  reflection  (R)  and  transmission  (T)  coefficients  are  given  by 
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R  =  b 
o 


T  =  ao  ' 


(45; 


Including  the  incident  TEM  mode  as  well,  the  dipole  model  of  the  scattering  obstacle  implicitly 
yields  a  symmetric  scattering  matrix  S  given  by 


(,:T  t) 


146! 


The   scattering   matrix   S   may   be   related   to   a   normalized   impedance   matrix   via 
Z  =  (I  +  S)  (I  -  S)~  where  I  is  the  identity  matrix.   The  equivalent  T-network  circuit,  shown  in 
figure  12a,  follows  [34]  from  the  impedance  matrix  Z.  The  impedances  Za  and  Z^  are  related  to  R  and 
T,  and  thus  aQ  and  bQ,  according  to 

7=7  ,  (R  -  T)  ,  =  7  r  a°  "  °  1 

a    c  (2  -  R  +  T)     cl2  +  a  -  b  > 

0     0 

(47) 


1\       1  ,     -,   , 1_   .  1^   1 

i  +  b 
o    o 


7        =  7   f — -)     -   —    7   =  -Z  ( = +  -1  -  -  2 

"b    c  lR  +  T  2>       2  la  c  la  +  b    2>       2  a 


where  Zc  is  the  characteristic  impedance  of  the  section  containing  the  scattering  obstacle  (in  this 
case  the  RCTL  impedance) . 

If  aQ  and  bQ  are  small,  which  is  typically  the  case,  then  Za  and  Z^  are  approximately  given  by 

(a  -  b  ) 
Z  --Z    °    ° 


a    c    2 

(48) 
L     -  -2  '   1   ^ 


b     cva  +  b  • 
o    o 


As   the   obstacle   size   vanishes,   both   a    and   b    tend   toward   zero.     Therefore 
Z  ->  o  (short  circuit)  and  Z,  ■>  =  (open  circuit)  and  the  equivalent  circuit  reduces  to  a  simple 
transmission  line  as  expected. 

As  an  example,  consider  a  spherical  conducting  obstacle;  in  this  case   a  and  a  are  given  by 
(18)  and  a  substitution  into  (44)  along  with  the  field  components  (34)  yields 
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Z       2 

a  =  -i  4u—  k  rf-1  F  L 
o        no  o  yaJ       xy 


Z        2 
b  =  -i  12k  —  k  r  (-)  F  ^ 


[49! 


2 
where  the  dimension! ess  quantity  F   defined  by 


F  2  =  I  i   S1'n;  yS  -  Py)  sin  Mx  sin  Ma  J  (Mg) 
xy      L  ,   sinh  M  b  o  3 

J        |_  mo=l 


[Y   cosh_ 
*■  ,   sin 
mo=l 


hM  6"^   cos  Mx  sin  Ma  J0^Mg^ 


(50: 
2 


contains  the  field  dependence  due  to  the  obstacle  location.  If  the  sphere  is  both  electrically 
small  (k  r  <  1)  and  small  compared  to  the  primary  cell  dimension  "a"  (r/a  <  1),  then  (48)  and  (49) 
give: 

Z        2 

Z  «  -i  Z  4n  —  k  r  (-)  F  ~ 

a      c    ti   o  va '        xy 


Z        2   „ 

Zk  »  -i  Z  {I6it  —  k  r  (-1  F  ] 

b      c  l    ti   o  laJ   xy1 


(51) 
■1 


Thus  we  see  that  the  loading  due  to  a  spherical  conductor  consists  essentially  of  a  small  series 
capacitance  ( Za  term)  and  a  large  shunt  capacitance  (Z^  term).  In  general,  conducting  objects 
placed  in  a  TEM  cell  should  behave  in  this  fashion. 

The  equivalent  T-network  representing  the  obstacle  loading  may  be  included  in  the  overall  TEM 
cell  transmission  line  circuit  as  shown  in  figure  12(b).  The  loading  is  assumed  to  be  centered  in 
the  RCTL  section  of  length  2L^CTL.  Ideally,  the  tapered  sections  of  length  LTApER  should  have  the 
same  impedance  as  the  RCTL,  i.e.  Z^p^  =  Zc;  however,  since  in  practice  this  is  generally  not  the 
case,  the  taper  sections  are  shown  explicitly. 

Tests  were  performed  in  a  typical  NBS  cell  to  assess  the  validity  of  the  above  equivalent 
circuit  representation.  The  cell  has  dimensions  a  =  15  cm,  b/a  =  1,  g/a  =  .17  which  yields  a 
theoretical  RCTL  impedance  of  51.6  Q  (Zc).  The  taper  and  RCTL  impedances  were  also  measured  using  a 
time  domain  ref lectometer  (TDR).  This  technique  yielded  1  «  51  Q    and  Z  «  49.5  Q.  Agreement 

between  the  theoretical  impedance  computations,  measured  impedances,  and  the  desired  50  Q  value  is 
usually  considered  acceptable  if  within  an  agreement  of  2  Q.  In  evaluating  the  TEM  cell  circuit 
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model,  the  measured  impedances  will  be  used.  The  RCTL  section  has  a  length  2L.,rT|  =  30  cm;  however, 
the  taper  length  is  somewhat  ambiguous  since  it  may  be  measured  along  the  center  of  the  inner 
conductor  (L-t^per  =  15  cm),  along  the  outside  taper  edge  (Ly^pf_R  =  25  cm),  or  in  between.  Because 
this  cell  has  a  narrow  gap,  the  fields  should  be  strongest  near  the  gap,  therefore  we  will  take  the 
outside  dimension  as  the  proper  choice,  i.e.,  Ly^pjrp  =  25  cm. 

The  cell  was  terminated  with  a  50  Q  load  and  a  spectrum  analyzer  was  used  to  measure  the  input 
impedance.  Three  spherical  conductors  were  used  to  load  the  cell.  Thus  equation  (51)  will  be  used 
to  determine  the  loaded  cell  equivalent  T-network.  We  begin  first  with  an  empty  cell.  Figure  13 
shows  the  magnitude  of  the  input  impedance  for  both  measured  data  and  theoretical  values  based  on 
the  circuit  of  figure  12(b)  with  one  end  terminated  with  a  50  Q  load.  At  very  low  frequencies 
(<  100  MHz),  the  well  defined  load  impedance  dominates.  As  the  frequency  is  increased,  the  taper 
impedance  (51  Q)  and  the  RCTL  impedance  (49.5  Q)  become  important.  The  theoretical  curve  appears  to 
be  first  dominated  by  the  taper  (around  250  MHz)  and  then  by  the  RCTL  (around  450  MHz).  The 
measured  data  shows  a  similar  behavior  up  to  400  MHz  whereupon  the  results  start  to  deviate 
significantly.  A  limiting  factor  on  the  circuit  model  is  the  stability  of  the  various  impedances 
with  frequency.  Although  measured  impedance  values  were  used,  TDR  results  represent  time  averaged 
values,  in  this  case  over  a  spectrum  reaching  into  the  12  GHz  range.  In  fact,  at  specific 
frequencies  the  actual  RCTL  and  taper  impedances  may  deviate  up  to  1  Q  or  so  from  these  average 
values.  In  practice  this  means  a  VSWR  of  less  than  1.05,  which  is  acceptable  from  the  cell  usage 
point  of  view,  but  it  is  not  reasonable  to  expect  the  circuit  model  to  give  better  than  the 
qualitative  agreement  found  in  figure  13.  What  our  analysis  does  give  is  the  basic  variation  in 
impedance  magnitude  with  frequency. 

Figure  14  shows  the  magnitude  of  the  input  impedance  for  the  same  parameters  as  in  the  previous 
unloaded  cell  case,  except  that  a  conducting  sphere  of  4.35  cm  in  diameter  has  been  introduced.  The 
sphere  is  centered  in  the  test  chamber,  i.e.  x  =  0  and  y/b  =  .5  and  this  particular  sphere  occupies 
slightly  less  than  a  third  (29%)  of  the  chamber  height.  The  curves  show  little  variation  from  the 
unloaded  cell  case,  except  above  400  MHz  where  the  measured  data  show  a  decrease  in  the  input 
impedance.  Next  a  6.5  cm  sphere  (approximately  43%  of  the  chamber  height)  was  tested  and  the 
results  are  shown  in  figure  15.  Again  the  theoretical  curve  is  little  perturbed  from  that  of  the 
empty  cell  while  the  measured  data  show  a  further  decrease  in  input  impedance  for  the  upper 
frequency  range  (>  400  MHz).  Finally,  we  introduced  a  sphere  occupying  two-thirds  the  chamber 
height  (diameter  =  10  cm  or  67%  of  the  chamber  height).  The  cell  is  significantly  perturbed  as 
shown  by  the  results  in  figure  16.  In  this  case  the  actual  perturbation  (measured)  is  greater  than 
theoretically  predicted,  most  likely  because  of  significant  coupling  between  the  sphere  and  the  cell 
walls.  In  addition,  this  sphere  excites  the  cavity  resonance  associated  with  the  TEq^  mode  at 
around  400  MHz,  and  the  transmission  line  model  is  no  longer  a  reasonable  description  of  the  cell 
behavi  or. 

Together  these  curves  demonstrate  that  the  small  obstacle  theory  gives  a  good  initial  model  of 
the  loaded  TEM  cell  test  circuit.  Up  to  the  first  cell  resonance,  the  theory  reasonably  well 
predicts  both  the  frequency  dependence  and  the  variations  in  the  input  impedance.  Thus  the  model 
should  be  useful  in  estimating  the  VSWR  change  expected  when  loading  a  cell.  More  importantly,  as 
the  loading  does  not  cause  significant  variation  for  the  initial  two  spheres  (<  43%  of  the  test 
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chamber  height),  the  present  analysis  supports  the  common  practice  of  ignoring  the  loading  effect 
when  the  size  of  the  test  object  is  no  greater  than  one  third  of  the  chamber  height. 

A  second  concern  is  the  change  in  the  field  distribution  due  to  a  test  object's  scattering  of 
the  TEM  mode  and  its  excitation  of  higher-order  (evanescent)  modes.  Only  the  RCTL  section  will  be 
considered  and  again,  for  simplicity,  a  centrally  located  (x  =  0,  y/b  =  .5),  spherical  conducting 
object,  of  diameter  .4b  will  be  used.  The  modes  will  be  scattered  according  to  equation  (1)  with 
the  expansion  coefficients  given  by  (10).  The  incident  mode  is  assumed  to  be  a  forward  travelling 
TEM  mode,  thus  the  TEM  mode  expansion  coefficients  are  given  by  (49)  while  those  for  the  TE-^q  and 
TE-q  modes  are  determined  by  (10)  with  the  proper  field  expressions,  (A-18)  and  (A-28)  respectively, 
inserted.  The  TEg-^  mode  should  not  be  excited  if  the  diameter  of  the  sphere  is  less  than  half  the 
test  chamber  height.  Modes  beyond  these  first  three  are  not  expected  to  contribute  at  normal  TEM 
cell  frequencies. 

Rather  than  consider  the  field  distribution  at  any  specific  point,  we  will  investigate  the 
perturbation  to  the  TEM  mode  field  components  (E  ,  Hx)  at  the  test  chamber  center  line  (x  =  0, 
y  =  l/2b)  as  a  function  of  z,  the  longitudinal  distance  from  the  scattering  sphere.  These  two 
components  are  of  primary  interest  since  they  form  the  basic  test  field.  The  wave  impedance  E  /Hx 
will  also  be  considered.  Each  parameter  will  be  normalized  by  its  empty  cell  value  at  that  point, 
thus  a  result  of  unity  represents  a  zero  perturbation  due  to  the  scattering  sphere.  Figure  17  shows 
the  numerical  results  of  the  above  analysis  with  the  frequency  chosen  to  give  an  electrical  cell 
width  of  k  a  =  1.  As  may  be  seen,  there  is  no  perturbation  in  the  forward  direction,  however,  the 
greater  scattering  of  the  TEM  mode  in  the  backward  direction  fib  I  >  la  I)  combined  with  the  phase 
beating  between  the  incident  and  reflected  modes  yields  a  noticeable  change  for  z  <  0.  The  electric 
field  is  reduced  while  the  magnetic  field  is  enchanced,  thus,  the  wave  impedance  shows  the  greatest 
deviation.  If  the  contributions  from  the  individual  modes  are  examined,  it  is  found  that  the 
scattered  TEM  mode  essentially  accounts  for  all  of  the  perturbation  shown  in  figure  17.  If  the 
frequency  is  increased,  so  that  k  a  =  1.4  which  nears  the  TE-^q  cutoff  value  of  n/2,  then  as 
expected  the  perturbation  is  increased  as  shown  in  figure  18.  The  reflected  fields  are  greatly 
altered,  resulting  in  a  wave  impedance  change  of  greater  than  10%.  In  addition,  we  now  see  some 
variation  in  the  forward  fields,  although  the  beating  effect  is  significantly  less.  Again,  the 
variation  is  largely  due  to  the  scattered  TEM  mode.  These  figures  indicate  that  it  may  not  be  safe 
to  assume  that  the  fields  are  unchanged  in  the  presence  of  a  test  object.  Noting  that  the  6.5  cm 
scattering  sphere  in  figure  15  is  nearly  of  the  same  diameter  (normalized  to  b)  as  in  figures  17  and 
18,  it  appears  that  the  input  impedance  is  much  less  sensitive  to  loading  than  is  the  field 
perturbation,  especially  for  z  <  0. 

Taken  together,  the  above  results  indicate  that  it  is  reasonable  to  assume  that  the  cell's 
impedance  characteristics  are  unchanged  by  a  small  test  object.  However,  more  care  must  be  taken  if 
the  same  assumption  is  to  be  applied  to  the  field  distribution. 

6.  Conclusions 

This  report  has  dealt  with  the  question  of  how  an  electrically  small  discontinuity  changes  the 
characteristics  of  a  TEM  cell.  This  quite  simple  formulation  was  then  applied  to  two  specific 
examples;  (1)  the  dual  TEM  cell  and  its  coupling  characteristics,  and  (2)  the  loading  effect  of  a 
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test  object  in  a  TEM  cell.  In  each  case,  experimental  data  were  compared  to  theoretical  results  in 
order  to  verify  the  usefulness  of  this  particular  analytical  model.   It  was  found  that  the  dual  TEM 
cell  behaves  as  predicted  if  the  aperture  is  empty;  while  for  a  loaded  aperture,  results  indicate 
qualitative  agreement  as  well  as  which  coupling  parameters  (primarily  the  contact  resistance)  need 
to  be  controlled  if  better  correlation  is  to  be  achieved.  The  analysis  of  TEM  cell  loading  suggests 
that  the  fields  inside  the  cell  are  more  sensitive  to  the  scattering  effect  of  the  test  object  than 
is  the  cell's  transmission  line  behavior,  such  as  its  characteristic  impedance.  Roth,  however, 
should  not  be  significantly  perturbed  in  the  course  of  normal  TEM  cell  usage. 

In  addition  to  these  two  present  applications  of  small  aperture  (obstacle)  theory  to  the  TEM 
cell,  others  may  arise  in  the  future.  For  example,  it  has  been  suggested  [72]  that  a  small  TEM 
cell,  with  an  aperture  in  one  of  the  walls  parallel  to  the  septum,  could  be  used  as  a  low  frequency 
field  probe.   If  the  incident  field  could  be  specified,  the  above  theory  would  again  predict  how  the 
TEM  mode  would  be  excited  in  the  cell,  and  consequently,  how  much  power  the  cell  would  receive.  The 
advantage  of  such  a  TEM  cell  probe  is  that  the  fields  couple  directly  into  a  known,  shielded 
transmission  line  thus  circumventing  the  problem  matching  the  probe  antenna  and  the  inevitable 
coaxial  feed  cable.  Also,  because  aperture  coupling  is  a  very  common  EMC  concern,  a  standard 
aperture  coupling  probe  could  be  used  to  predict  the  impact  the  same  aperture  would  have  if  present 
on  a  piece  of  electronic  equipment.  Thus,  TEM  cell  probe  results  could  be  used  as  a  building  block 
toward  developing  more  complicated  EMC/EMI  models  for  electronic  systems. 
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Figure  2.     Dual    TEM  cells  with   common   aperture, 
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Figure  3.  The  cross  section  of  an  RCTL 
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Figure  4.   Received  power  in  the  sensing  cell  for 

the  case  of  well  joined  cells  about  a  square  aperture 
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Figure  9.      Insertion   losses  measured  at   the  forward  port    (P3)    for  a   layered  polyester-aluminum-vinyl 
sample    (after  thickness   correction)   versus   frequency,   compared  with  contact   resistance   (In- 
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38 


Hi 


5Q 
44 
33 
32 
25 
20 
H 

3 

2 

-4   h 


ME5H  PRTTERN  IL  -  BRCKWRRD  PORT  -  THICKNESS  CORRECTED 


-10 


T 


T 


T 


T 


T 


X 


.*■- 


*:--* 


RC  =  D 


RC-1 


RC-Z 


****    MEREURED    VRLUEB 


8  100         2@Q         300         408         500         E0D         780         80S         3  80         1088 

FREQUENCY    ( MHz ) 


Figure  10.    Insertion   losses  measured  at  the  backward  port    (P2)   for  a  layered  polyester-alumi num- 
vinyl    sample   (after  thickness   correction)   versus  frequency,  compared  with   contact   resistance 
(Rc)-dependent  theoretical    curves. 
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Figure  11.    Insertion   losses   due  to  a   layered  polyester-aluminum-vinyl    sample  measured  with   various 
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(a)  The  equivalent  T-network  circuit  for  a 
small  scattering  obstacle. 
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(b)  The  TEM  cell  transmission-line  circuit, 
including  the  small  obstacle  loading. 


Figure  12.  TEM  cell  equivalent  circuit  representations 
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Figure  13.  Input  impedance  versus  frequency  for  the  empty  TEM  cell 
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Figure  14.  Input  impedance  versus  frequency  for  the  TEM  cell  with  a  4.35  cm  diameter  conducting 
sphere  present. 
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Figure  15.    Input   impedance  versus   frequency   for  the  TEM  cell   with   a  6.50  cm  diameter  conducting 
sphere  present. 
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Figure  16.  Input  impedance  versus  frequency  for  the  TEM  cell  with  a  10.00  cm  diameter  conducting 
sphere  present. 
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Figure  18.  Field  perturbation  due  to  a  conducting  sphere  of  diameter  6.0  cm  with  kQa  -  1.4, 
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APPENDIX  A 

The  Normalized  Modes  In  a  TEM  Cell 

The  higher-order  modes  in  a  TEM  cell  are  difficult  to  analyze  and  are  largely  ignored. 
Primarily,  interest  is  in  the  dominant  TEM  mode  and  the  associated  RCTL  impedance  characteristics. 
The  higher-order  modes  are  considered  only  to  the  extent  that  their  appearance  limits  the  TEM  cell's 
usable  bandwidth.  Thus,  given  the  current  TEM  cell  applications,  a  complete  modal  description  is 
both  unavailable  and  unnecessary.  The  present  investigation  will  only  look  at  the  first  four 
hi  gher-order  modes . 

Several  solutions  for  the  TEM  mode  fields  have  appeared.  For  the  case  of  a  centrally  located 
inner  conductor,  a  conformal  mapping  approach  [21]  accurately  describes  the  fields  in  terms  of 
elliptic  functions.  Alternately,  the  singular  integral  equation  method  [23],  along  with  the 
condition  that  the  gaps  (see  fig.  3)  be  small,  results  in  field  expressions  involving  infinite  sums 
of  simple  functions.  The  latter  solution  allows  for  a  vertically  offset  inner  conductor  and  also 
covers  the  higher-order  modes  of  interest  here.  Therefore,  the  singular  integral  equation  results 
will  be  used  here,  although  strictly  speaking,  this  solution  is  not  accurate  for  points  very  near 
the  gaps.  Accurate  results  are  obtained  in  the  central  test  area,  however. 

The  first  four  RCTL  modes  are  the  TEM,  TEqj,  TE-^q,  and  TE-jj  modes  respectively  [73]  for  typical 
cell  dimensions  (small  gap).  However,  the  TEq^  mode  does  not  interact  with  a  test  object  centrally 
located  in  the  test  zone  (x  =  0)  as  this  mode  largely  confined  to  the  gaps  [74].  Therefore  the 
TEqj  mode  may  be  ignored.  The  other  three  modes  need  to  be  orthonormal ized  according  to  equation 
(4).  For  m  *  n,  (4)  is  satisfied  via  modal  orthogonality  (34);  while  for  m  =  n,  (4)  reduces  to 

/  (E2  +  E2  )  ds  =  Z  (A-l) 

rj.  *■  nx    ny J      n  v 

where   Zn  is  the  wave   impedance  of  the   n-th  mode. 

The  desired  TE-modes  may  be  found  from  a  knowledge  of  each  mode's  longitudinal  magnetic 
field.  The  singular  integral  equation  approach  mentioned  above  solves  for  the  TE-fields  excited  by 
an  elementary   vertical    electric  dipole.      In  the  transform  domain    (z-di rected) ,   we  find  that 


cos   K   (b  .   -  py ) 

mv  j       VJ  J 

K     s  i  n  K     b  . 

mo         m  m     j 


HU)=p^-a      (p)      I    ^ —  sin  Mx   sin  Ma  J     (Mg) 

z  K  a       o     H       L         KsinKb.  o9 


Id!  G.{x,   y,   0,   yQ)    ,  (A-2; 


where    the    RCTL    dimensions    are    defined    in    figure    3,    p    is    the    transform    variable    corresponding   to 

propagation       in       the       z-direction,       j       =       1,2       refers       to      the      upper      and       lower      chambers 

7  7  7 

respectively,  Km  =   (M     -  k~j  +  p   j        ,    Id!  is  the  dipole  moment,    (0,  y0)   is   the  dipole  location,   and 
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t^o  o,  -1  M  cos   K      (b,    -  y   ) 

33  mo  m  m     1 

.2  cot   K  b .       , 

j=l       mo  m 

,  M  cos   K      (b .    -  py) 

Gj    (X'  y'    °»   V    =  HT       Z  K     sin   KJ   b.        Sin   Mx   cos   Km  yo 

mo  m  m     j 

The  remaining  notation  is  defined  in  conjunction  with  (34).  Note  that  the  dipole  excitation  nature 
of  these  equations  is  extraneous  here,  however,  for  simplicity  the  results  of  [23]  will  be  applied 
di  rectly . 

The  TE  modes  may  now  be  found  by  evaluating  the  appropriate  pole.  We  may  first 
evaluate  H  l  which  will  determine  the  remaining  TE  field  components,  or  we  may  use  the  following 
relationships 


~    (j)=  ionx  ~   (j)  H~(JL^±ora     H~   (J) 

k2   _  a2     y     z  x  2         2     x     z 

O  0 

(A-4) 

E    (j)=        9i("%    5     H   (j)  H     =        :ia     9    5     H   (j) 

y  k   2   -   a2     X      Z  y        k   2    -   a2     y      Z 

0  0 

and  again  evaluate  the  proper  residue.  One  final  comment  about  H  J  is  in  order  before  proceeding 
to  specific  modes.  The  first  term  in  (A-2)  contains  the  gap  perturbation  information  and  therefore 
is  associated  with  gap  dependent  modes  such  as  the  TEM  mode.  If  there  is  no  gap,  the  TEM  mode 
cannot  exist.  The  second  term  is  an  unperturbed  term  associated  with  the  modes  which  would  exist  in 
the  dipole  chamber  if  the  gap  were  closed.  Thus  this  term  will  not  contribute  to  the  TEM  or  TEii 
modes  . 

(1 )   TEM  mode:  the  TEM  mode  is  due  to  poles  at  p  =  ±  k  .As  mentioned,  the  second  term  in  (A-2) 

°  (i) 

does  not  contribute  (further  details  may  be  found  in  [23],  Appendix  C).  The  H  J  term  will  be  zero 

for  the  TEM  mode,  thus  the  field  components  need  to  be  found  directly  from  (A-4).   Because  we  are 

after  the  field  distribution,  we  need  only  evaluate  one  pole,  say  the  forward  propagating 

mode  (p  =  -k  ),  and  consider  the  electric  field  components  since  the  magnetic  field  components  are 

related  via  the  intrinsic  wave  impedance  ti  .  A  residue  calculation  yields 

t)  ng  sinh  M(b  .  -  py) 

E  =  ~ a  (-k  )   I    r-u-Eik r  sin  Mx  sin  Ma  J   (Mg) 

x    2a    oK     o    L         sinh  Mb.  o   3 


mo 


p  ri  -ng  cosh  M(b  .  -  py) 

E  =  -~ a  (-k  )  I    .  .  I. cos  Mx  sin  Ma  J  (Mg! 

y     2a     ov  <r  L  sinh  Mb.  o   3' 

mo         j 
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:a-5: 


From  [23,  eq .  58]  we  find  that 


a„(-kJ  =  4^—  (A-6) 


where  V  is  the  TEM  mode  voltage.  Thus  (A-5)  may  be  written 


sinh  M(b  -  py) 

E  =  t±     )'  .  JM, sin  Mx  sin  Ma  J  (Mg) 

x   a   L  sinh  Mb.  o   3 

mo  j 


2y  cosh  M(b.  -  py) 

E  =  p —  )  .  r J... cos  Mx  sin  Ma  J  (Mg) 

y   'a  L  sinh  Mb.  o   a 

J  mo  j 


:a-7) 


To  normalize  we  note  that  the  power  (P)  carried  by  the  TEM  mode  is  given  by 


2 
/  (E  ±  x  H  **)  •  a"  =  +  2P  =  +  ¥=-   .  (A-8) 

CS    °     °      z  Zc 


If  we  expand  the  integral    via    (2)   we  find  that  (A-8)   is  equivalent  to 

/  TEot  x"Hot)    .az-£     .                                                            (A-9) 
CS                                               c 

1/2 

Therefore,  we  need  only  normalize  by  Z  /V  to  satisfy  (4)  or  alternately  (A-l).  Thus  we  arrive  at 

y       1/9  sinh  M(b.  -  py) 

E  =  -  I   L,c      I    ■  r  J sin  Mx  sin  Ma  J  (Mg) 


"x   a  c     i    sinh  Mb .         ..-     , 
mo 


(A-10) 


?   1/?    cosh  M(b,  -  py) 

E  =  p^  Z  '        ]    .  .  J cos  Mx  sin  Ma  J  (Mg) 

y   Ka  c     ^    sinh  Mb .  o  v  3 

J  mo        j 

The  magnetic  components  are  determined  by  (9). 

(2)    TE^q  mode:  the  TE-^g  mode  is  due  to  poles  at  p  =  +  plf)  in  both  the  perturbed  and 
unperturbed  terms  of  (A-2),  where  p,n  is  the  unperturbed  rectangular  waveguide  TE^q  mode  propagation 

constant,  i.e.   p1Q  =  (kQ2  -  (rc/2a)2)  1/2 .   Note  that  for  p  =  +  p1Q,  Kj  =  0.  For  K1  small,  (A-2) 

yields 
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p?ao<fW 


si  n 


7IX 

2a 


b.    K, 
J      1 


2     o    va;  ij 


KX 

Id  A  n    Sin  2a 
a     2a  bx   K}2 


(A-ll) 


ao  {ho' 


■Id  A     2      ,      >iig> 
g     4ab     o    l2aj 


where   we    have    let    the   dipole    location    (y   )    be    in    the    upper    chamber    (j=l).      The    proper    residue    for 
the  forward   TE-^q  mode  follows   readily   from   (A-ll)   and  we  find 


■   yj)   _  -.-Ida     it     drl   nx    ,6ij         P  2      ,   2   ,jg_n 

■  ]7T  6^sin  2?  {"¥T-  2FT7  Jo     (2T)}  * 
1  J 


4a2   P10 


[A-12] 


Notice    that    for    j  =  l,2    we    find    that    the    term    in    parenthesis     gives     l/2b    if    terms    on    the    order 
of   (iig/2a)    .   are  neglected.     Thus 


u    (J)    -   j  IdJL       ti       .      itx  _    .      ,      nx 

H  -   i — y-  - —  sin  t-  -  A  sin   j- 

z  8a2b   ho  2a  2a 


[A-13] 


As  expected,  (A-13)  simply  describes  a  rectangular  waveguide  TE-^q  mode  with  the  magnitude  A 
determined  by  the  dipole  strength  and  guide  dimensions.  The  electric  field  components  follow  from 
(A-13),   i.e. 


"X 
:(Ji 


0 


.,  2a   .  tix 

ik   -n    —  A  cos  7j- 

o'o  1  2a 


:a-h) 


We  now  apply  the  normalization  condition  (A-l)  which  requires  that 


(■    1    2a^2  .2  t       „  2  nx    . 

h  k  ti  —   A   J  cos  -7T-  ds  =  Z 

1   0  'o  it  J     r^      2a 


10 


(A-15) 


where  the  TE-^q  wave  impedance  is  given  by, 


Z10  =  koVp10  * 


:a-i6] 


Combining     (A-15)     and    (A-16),     integrating    over    the    cross    section       (x    e  [-a,    a],  y    e  [-b~,   b,]),  we 
find   that 


zin  1/2 


2ab;         i    k   ti     2a 
0  0 


[A-17) 
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Thus  the  proper  normalized  components  are 


EVJ/  =  0 
x 

E(j)  =  (  *°)    cos^ 

y  vab;       2a 

/  ■>      •     Z1n  1/2 

(J)  ■  -^-(oS    sin- 


IA-18] 


H 


2a  k  n  ^2ab 
o  'o 


2a 


The  remaining  magnetic  components  follow  from  (9).  The  final  result  (A-18)  could  have  been  written 
directly,  since  for  the  TE^q  mode,  the  RCTL  is  nothing  but  a  rectangular  waveguide  of  cross  section 
2a  by  2b.  However,  developing  (A-18)  via  the  full  residue  calculation  serves  as  a  useful  guide  when 
evaluating  the  more  complicated  TEii  mode.  In  the  case  of  a  zero  gap,  the  TE^  mode  should  reduce 
to  the  TE1Q  result  (A-18)  . 

(3)   TE-q  mode:   the  TE^  mode  arises  from  poles  in  a  ( p)  when   p  =  +  p., ,  ,  where  p, ,  ->  p,  „  as 
g  +   o.   The  unperturbed  term  will  not  contribute.  Again  evaluating  the  forward  mode  (p  =  -Pii)  we 
find 


, . >  cos   K     (b  .   -  py) 

HU;   =  pMResidue(a    (p     ))      I  m. — I — r sin  Mx  sin  Ma  J      (Mg)    . 

z     Ka         o  Hll    L         K  sin  Kb.  o   3 

mo    m     m  j 

(A-19) 

Because  we  expect  the  TE-q  mode  to  propagate  nearly  the  same  as  the  TE^q  mode,  we  may  write  p, ,  as 
fol lows 


2     2     2 

311  =  P10  "  611 


'A-20) 


where  6,,   should  be  small 


2  2 

6,,    ■>  0  as   g  ->  0).   Indeed,   if     6;,    <<   1,  then 


6n  m=l 


iM  m>3 


:A-21) 


Substituting   (A-21)    into   (A-19)  yields 


i(j)   -  ~^ 


cos   6n(bj   -  py! 


HlJj   -  p2l  Residue(a    ( p"  .  ) )  ,,,,   ,.     „      ,  , 

oK*  i  i         I     6,,    sin   6,,    b.  2a     o    v 2a 


tcx    ,      /  irg  \ 
sin  —J     (of) 


=°  cosh   M   (b  .   -  py) 

y        53 — ■   ,    M, sin  Mx   sin  Ma  J      (Mg' 

L  -,  M  sinh    Mb .                                         o    v    b' 

mo=3  j 


(A-22) 
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To  normalize  the  TE^  mode  we  may  ignore  the  terms  outside  the  brackets  as  was  done  with  the 
TE1Q  mode  in  (A-13).  However,  to  insure  that  hTJ'  tends  toward  the  proper  limit  as  g  ■*■  0  an 
additional    term  of 


J0Ug/2a)    /(6n   sin   &u    b . ) 


will    be  factored  out.     Thus  the   relevant   components   are  given  by 
,(j)   _ 


H' 


p   A   {cos   6, , (b  .   -  py)    si  n 


T[X 

2a 


sin   6,,    b.        cd       cosh   M(b  .   -  py) 

-6.,  -=— i       I  .JM, sin   Mx  sin  Ma  J      (Mg)    } 

11      ,      /itgi  L  -  sinhMb.  oa 

J0    {■£)  mo=3  j 


Exj)   =   "iko%   &     611  A   tsin   6n(bj   -py)   sin-g 


si  n  6, ,    b  .        <*>       sinh   M(b  .    -  py) 

+ ±i_J.      y __J sin 


/■ng->  ^  _'       sinh  Mb . 

Jo   (2?)       mo=3  J 

2a  2 

1  k0%  (~}  6n  p  A  {2f\[  cos  6n(bj  -  py)  cos 


Mx  sin  Ma  J  (Mg) } 


■nx 
2a 


:A-23) 


si  n  6, 1  b  .   co  cosh  M(b  .  -  py) 

— ~       I  ■  c  mk cos  Mx  sin  Ma  J  (Mg)} 

,  r iig a     L   o    sinh  Mb .  0   a 

J0  (2?)  mo=3         j 


Note  that  as  6,,  ■*■  0  (no  gap),  the  field  components  reduce  to  the  TEjq  [expressions  (A-13)  and  (A- 
14)]  with  only  an  additional  p  factor  which  accounts  for  the  even  (TE^q)  or  odd  (TE-q)  polarity  of 
the  TE-^q  modes  in  the  individual  chambers. 

These  expressions  may  be  substituted  into  (A-l).   Noting  the  orthogonality  between  cos  Mx  and 
sin  Mx  over  the  interval  (-a,  a),  one  finds  the  following 


4  1 

*2  ,2a^   c2  „2 


py; 


(1  ko%"  O  6ii A  a  _/     fe2^  cos  an<bj  -  ^  +  sin  6n(bj 


sin  6,,  b .  co    j  (Mg) 

-2     "  J   I       -2-2 cosh  2M(b  -  py)  dy  =  Z 

iT  f-29-)  mo=3  sirrtr  Mb.        J 

0  v2a;  .-1 


where  Z-q  is  defined  as  in  (A-16) .  Evaluating  the  remaining  integral  yields 


'A-24; 
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:'  koV2  'f'4  6n  A?  a     I    (<?ri^»2  d^  (611  bj  +  1  s1"  611  b; 


jj-  (6U   bj    -j  sin   2   6n    bj)   ♦ 


2  2 
sin     6, ,    b.  =o      J        (Mg) 
11     J        y         o     v    a 

i   2 ,  -rig  ^  ^  _         M 

JQ    (-£)  mo=3 


coth   Mb.  }  =   Z,. 


IA-25! 


4  2 

If  we  now  neglect  terms   like   6,,,  and    (-n:g/2a)",  then    (A-25)    reduces  to 


!1    k0.o)2    (f-f   A2   2ab   =   Zn    . 


(A-26! 


Thus  we  find  that 


1/2 


11 


l2abJ    i  k  n  2a 
o  o 


(A-27) 


which     is    equivalent,    as    might    have    been    expected,     to     (A-17). 
normalized  components: 


Combining     results,    we    find    the 


(j)     .         (Hi  r  11     _, 


baB*        {_T—  sin  5ll  (bj  -  py)  sin 


TIX 

2a 


2a  6,, 

+  —  sin  6,,    b.        I 

tc  11     J     ^ 


»       sinh   M(b  .   -  py) 

.   ,    .,1 sin   Mx  sin  Ma  J      (Mg)} 

sinh  Mb .  o   v    3 


mo=3 


Z       1/2 
Eyj)   -^  {C0S   6ll(bj   "P^    C0Sf 


:a-28] 


2a  6 


11 


sin  6, ,   b  .       V 

71  11     J     _£ 


cosh   M(b  .   -  py] 


-  cos   Mx  sin  Ma  J       Mg) } 
;         u  0  sinh   Mb .  0    v    3/ 

'     mo=3  j 


^J)   -         P 

z  i   k     n      ^2a 

o    'o 


Z       ^2 

)    ^2ab^  {cos   6ll(bi   "  py^   sin  I 


cosh  M(b  .  -  py' 

I 
mo=3  j 


>, ,  sin  6,,  b.   V   — u — ■  .  \,.         -sin  Mx  sin  Ma  J  (Mg)} 
11      11  J    _o   M  sinh  Mb,  o  v 


where  the  small  gap  approximation  has  been  retained.   Again,  note  that  as   6,,  *  0  (no  gap)  these 
reduce  to  the  proper  TE^y  forms  modulo  polarity. 

Finally,  for  the  small  gap  case  considered  here,  it  is  found  [16],  [23]  that   6,,   is 
approximately  given  by 
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6, ,    =  it  b,  b„  /ab   Zn[ — 1 
11  12  v-nn; 


:A-291 


Tnis  result  may  be  used  to  determine  p, ,  and  Z^j  via  (A-20)  and  (A-16)  respectively.  Actually, 
(A-29)  is  not  a  particularly  accurate  expression  since  the  assumption  of  a  small  gap  implies 
that  In  [8a/ng)  >>  1,  which  is  not  the  case  even  for  g/a  =  .1.  A  better  approximation  follows  from 
solving  the  TE-,i    small    gap   modal    equation    [16,    23] 


o — j—-  (cot   6, ,    b,    +  cot    6,,    b„)    =    Jin  f — 1  -2 
2a    6,  ,  111  112  '"nq ; 


11 


;a-3o) 


from  which  (A-29)  was  found  assuming  both  that  6,,  b.  <<  1  and  Jin  (8a/itg)  >>  1 
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teristics of  building  elements  and  systems. 

Technical  Notes — Studies  or  reports  which  are  complete  in  them- 
selves but  restrictive  in  their  treatment  of  a  subject.  Analogous  to 
monographs  but  not  so  comprehensive  in  scope  or  definitive  in 
treatment  of  the  subject  area.  Often  serve  as  a  vehicle  for  final 
reports  of  work  performed  at  NBS  under  the  sponsorship  of  other 
government  agencies. 

Voluntary  Product  Standards — Developed  under  procedures 
published  by  the  Department  of  Commerce  in  Part  10,  Title  15,  of 
the  Code  of  Federal  Regulations.  The  standards  establish 
nationally  recognized  requirements  for  products,  and  provide  all 
concerned  interests  with  a  basis  for  common  understanding  of  the 
characteristics  of  the  products.  NBS  administers  this  program  as  a 
supplement  to  the  activities  of  the  private  sector  standardizing 
organizations. 

Consumer  Information  Series — Practical  information,  based  on 
N  BS  research  and  experience,  covering  areas  of  interest  to  the  con- 
sumer. Easily  understandable  language  and  illustrations  provide 
useful  background  knowledge  for  shopping  in  today's  tech- 
nological marketplace. 

Order  the  above  NBS  publications  from:  Superintendent  of  Docu- 
ments, Government  Printing  Office.  Washington.  DC  20402. 
Order  the  following  NBS  publications— FI PS  and  NBSlRs—from 
the  National  Technical  Information  Service,  Springfield,  VA  22161. 

Federal    Information    Processing    Standards    Publications   (FIPS 

PUB)— Publications  in  this  series  collectively  constitute  the 
Federal  Information  Processing  Standards  Register.  The  Register 
serves  as  the  official  source  of  information  in  the  Federal  Govern- 
ment regarding  standards  issued  by  NBS  pursuant  to  the  Federal 
Property  and  Administrative  Services  Act  of  1949  as  amended. 
Public  Law  89-306  (79  Slat.  1127),  and  as  implemented  by  Ex- 
ecutive Order  1 1 7 1 7  (38  FR  1 23 1 5,  dated  May  1 1 ,  1 973)  and  Part  6 
of  Title  15  CFR  (Code  of  Federal  Regulations). 

NBS  Interagency  Reports  (NBSIR)— A  special  series  of  interim  or 
final  reports  on  work  performed  by  NBS  for  outside  sponsors 
(both  government  and  non-government).  In  general,  initial  dis- 
tribution is  handled  by  the  sponsor;  public  distribution  is  by  the 
National  Technical  Information  Service  ,  Springfield,  VA  22161, 
in  paper  copy  or  microfiche  form. 
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